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Precision spectroscopy of simple atomic systems has refined our understand- 
ing of the fundamental laws of quantum physics. In particular, helium spec- 
troscopy has played a crucial role in describing two-electron interactions, de- 
termining the fine-structure constant and extracting the size of the helium nu- 
cleus. Here we present a measurement of the doubly-forbidden 1557-nanometer 
transition connecting the two metastable states of helium (the lowest energy 
triplet state 2 ^Si and first excited singlet state 2 ^Sq), for which quantum elec- 
trodynamic and nuclear size effects are very strong. This transition is fourteen 
orders of magnitude weaker than the most predominantly measured transi- 
tion in helium. Ultracold, sub-microkelvin, fermionic ^He and bosonic ^He 
atoms are used to obtain a precision of 8x10^^, providing a stringent test of 
two-electron quantum electrodynamic theory and of nuclear few-body theory. 
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The first observations of lielium emission spectra at the end of the 19th century revealed 
two separate series of lines, associated with orthohelium and parahelium respectively. In 1926 
Heisenberg explained the distinction between these two spectra on the basis of wave mechan- 
ics, electron spin and the Pauli exclusion principle (1). The spectrum of orthohelium arises 
from triplet states for which the electron spins are parallel, whereas in parahelium the electron 
spins are anti-parallel, forming singlet states (Fig. [T]). From the lowest state of orthohelium, 
the ls2s ^Si state (denoted 2 ^Si), only excitations to triplet states have been observed. Or- 
thohelium transitions from the 2 '^Si state and associated studies of the n ^Po,i,2 {n=2,3) fine- 
structure splittings (2-7) have enabled tests of quantum electrodynamics (QED) {8, 9) as well as 
a determination of the fine-structure constant (5, 10). In the singlet spectrum of helium (parahe- 
lium), electric-dipole transitions from the 1 ^Sq ground state {11) and from the metastable 2 ^Sq 
state {12, 13) have also provided tests of high-precision QED calculations. All these frequency 
metrology studies have been performed using either atomic beams or gas discharges. However, 
helium in the metastable 2 ^Si state (He*, lifetime 8 x 10'^ s) can be laser cooled and trapped, 
which allows much longer interaction times for excitation of weak transitions. He* atoms have 
been cooled to /iK temperatures, revealing quantum statistical effects of bunching and anti- 
bunching {14) and achieving quantum degeneracy for both the bosonic isotope ''He {15, 16) and 
the fermionic isotope ^He (77). 

Here we observe an orthohelium - parahelium transition, specifically, the 1557-nm transi- 
tion between the metastable 2 '^Si and 2 ^Sq states (Fig.[T]), both in ''He and '^He. This transition 
is an excellent testing ground for fundamental theory of atomic structure. Due to a large elec- 
tron density at the nucleus, the energy of S-states is the most sensitive to QED and to nuclear 
size effects {8). For the 2 ^Si and 2 ^Sq metastable states, QED terms contribute 4 and 3 GHz 
respectively, to a total binding energy of 10^ GHz {8, 9). The present accuracy in the QED cal- 
culations is 2 MHz, based on an estimate of non-evaluated higher-order terms. Many of these 
terms are common between the isotopes. Therefore, in the calculation of the isotope shift (i.e., 
the difference between the transition frequencies for ''He and '^He) mass-independent terms can- 
cel and the uncertainty is reduced to the sub-kHz level {18). As the finite nuclear charge radius 
shifts the 2 '^Si state by 2.6 MHz and the 2 ^Sq state by 2.0 MHz, an accurate isotope-shift 
measurement allows a sensitive determination of the difference in the mean charge radius of the 
a-particle and of the ^He nucleus, providing a stringent test of nuclear charge radius calculations 
and experiments {19). 

The natural linewidth of the 2 '^Si 2 ^Sq transition is 8 Hz, determined by the 20-ms 
lifetime of the 2 ^Sq state which relaxes via two-photon decay to the ground state (Fig.[T]). This 
transition is 200 000 times narrower than the natural linewidth of the 2 '^P state, which is most 
prominently used for spectroscopy in helium. The Einstein A-coefficient for the 2 '^Si— t- 2 ^Sq 
magnetic-dipole transition is ~10^^ s^^ {20, 21), fourteen orders of magnitude smaller than 
for the electric-dipole transition from 2 ^Si to 2 '^Po,i,2 states, which indicates that excitation 
requires high power and/or long interaction times. 

The experiment described here was performed using an apparatus designed for the pro- 
duction of quantum degenerate gases of helium {17, 22). Briefly, the metastable 2 ^Si state is 
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populated by electron impact in an electric discharge. The atomic beam is coUimated, slowed 
and trapped using standard laser cooling and trapping techniques on the 2 '^Si — )• 2 '^P2 transition 
at 1083 nm. The atoms, optically pumped to mj=+l, are then transferred to a loffe-Pritchard 
type magnetic trap. "^He* atoms are evaporatively cooled towards Bose-Einstein condensation 
by stimulating radio-frequency (RF) transitions to untrapped states. For "^He* (in the F=3/2 
hyperfine state), quantum degeneracy is reached by sympathetic cooling with "^He*. Either one, 
or both, of the two isotopes are transferred into a crossed-beam optical dipole trap. This trap 
consists of two focused 1557-nm laser beams, intersecting at their foci as shown in Fig.[2j We 
transfer up to 10^ atoms to this optical trap. 

After loading the optical trap, the atoms are illuminated by a separate beam for spectroscopy 
which is derived from the same laser as the optical trap beam, but is switched and frequency 
shifted by a 40-MHz acousto-optic modulator. A heterodyne signal is set up between the 1557- 
nm laser and a mode of a femtosecond frequency comb laser to deduce the absolute frequency 
of the spectroscopy laser. The frequency comb is based on a mode-locked erbium-doped fiber 
laser, for which both the repetition rate and the carrier-envelope-offset frequency are referenced 
to a GPS-controUed Rubidium clock (23). 

After a certain interaction time (typically 1 to 6 s), both the spectroscopy beam and the trap 
beam are switched off, allowing the atoms to fall due to gravity. The high internal energy of 
He* (20 eV above the 1 ^Sq ground state) allows for efficient detection on a microchannel plate 
(MCP) detector (Fig. |2]). The MCP signal reflects both the number of atoms and their tempera- 
ture. In the case of ^He, the signal has a bimodal character that results from the combination of 
Bose-condensed atoms and thermal atoms (Fig. 3(a) I; a fit to this signal provides the number of 
condensed atoms (23). Because the excited state is anti-trapped, the trap is depleted when the 
spectroscopy beam is resonant with the atomic transition. By deducing the remaining number 
of 2 ^Si atoms for various laser frequencies, the atomic resonance frequency is determined from 
a Gaussian fit to the data (Fig. |3(b)| ). The observed hnewidth is largely due to the 75-kHz laser 
linewidth. 

Several systematic shifts in the transition frequency are taken into account (23). The largest 
shift is due to the Zeeman effect. The measured transition, 2 '^Si (mj=+l) -^^ 2 ^Sq (mj=0) 
for ^He, and 2 ^Si, F=3/2 (mF=+3/2) 2 ^Sq, F=l/2 (m^=-i-l/2) for ^He, is shifted from 
resonance predominantly by the Earth's magnetic field. The size of the shift is deduced by 
measuring the resonance frequency of RF spin-flip transitions between the 2 ^Si magnetic sub- 
states. An additional shift is caused by the momentum transfer from a 1557-nm photon to an 
atom. In the case of ^He, the high density of the condensate could potentially cause a mean- 
field shift (24). However, by performing the experiment with reduced atomic density, no shift 
is observed. 

The second largest systematic frequency perturbation is due to the AC Stark shift associ- 
ated with the intense 1557-nm light which induces the dipole trap: the specific energy state of 
the trapping potential for an atom determines the AC Stark shift for that atom. For ^He, only 
excitations of atoms condensed in the ground state of the dipole trap are taken into account 
in determining the transition frequency. As the trap depth depends linearly on laser intensity. 
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measuring the resonance frequency for a range of applied laser powers allows an extrapolation 
to zero laser intensity. In contrast, ^He atoms, due to their fermionic nature, are distributed 
throughout the energy states of the dipole trap and as a result the measured AC Stark shift does 
not equal the trap depth (as is the case with ''He), but is reduced due to the density of states 
within the dipole trap. A non-linear shift can then potentially arise at high laser intensities, 
where the larger trap depths allow for higher temperatures. To minimize this effect "^He atoms 
are sympathetically cooled to the quantum degenerate regime to predominantly populate the 
lowest energy states of the trapping potential. Over the course of several months, 20 indepen- 
dent extrapolations were obtained (as shown in Fig. |4]) to deduce an absolute frequency of the 
2 ^Si ^ 2 ^So transition for ^He of /4 = 192 510 702 145.6(1.8) kHz and for ^He (F=3/2 
-> F=l/2) of /a = 192 504 914 426.4(1.5) kHz, where the one standard deviation error in 
parentheses includes all statistical and systematic uncertainties. 

For both isotopes, our result agrees with QED calculations of the ionization energies of 
the two metastable states (9, 25). The present experimental error in the transition frequency is 
three orders of magnitude smaller than estimates of non-evaluated higher-order terms in state- 
of-the-art QED calculations and presents a significant challenge for groups involved in atomic 
structure theory. 

An indirect value of the energy difference between the 2 ^Si and the 2 ^Sq states can be ob- 
tained from the literature (only for '^He) by combining experimental transition frequencies from 
both metastable states to high- lying S, P and D states with theoretical values for the ionization 
energies of these states. This procedure yields ionization energies for the 2 ^Sq state (8, 12, 13) 
and the 2 ^Si state (2, 8) and the difference between these values gives a transition frequency of 
192 510 701.96(16) MHz, in agreement with our result though with hundred-fold lower preci- 
sion. 

Isotope shift measurements, combined with high-precision QED theory, provide a method 
to isolate contributions due to finite nuclear size effects. The difference in nuclear charge radii 
between ^He and ^He is determined by comparing experiment and theory. Because the ^He 
nuclear charge radius is the most precisely known radius of all nuclei determined from elec- 
tron scattering experiments (2(5), 1.681(4) fm, a value of the '^He nuclear charge radius with 
similar precision can be deduced. In calculating the isotope shift, QED theory is more precise 
than our measurement as mass -independent terms cancel. The theoretical value for the iso- 
tope shift (assuming point-like nuclei) is 8 034 148.6(7) kHz (23). Subtracting the measured 
transition frequencies and correcting for the accurately-known hyperfine structure (25, 27), we 
find an isotope shift of - h + Ufs = 8 034 367.2(2.3) kHz. The 218.6 kHz difference 
may be attributed to the finite size of both nuclei. This nuclear shift is proportional to the 
difference in the nuclear charge radii squared, Ar^ = ("^He) — (^He). Using the theoret- 
ical proportionality constant of 4.6642 fm^/MHz (18) for the measured transition, we deduce 
Ar^ = 1.019(11) fm^. Ar^ represents a more universal parameter than the value of the isotope 
shift as it is obtained from various branches of physics. Besides through spectroscopic means, 
it can be determined from nuclear theory and from electron-scattering experiments. Nuclear 
few -body theory provides Ar^ = 1.16 ± 0.12 fm^ (18,23,28), whereas from electron-scattering 
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Figure 1 : Relevant energy levels, transition wavelengths and state lifetimes of helium. The 
magnetic-dipole transition connecting the 2 ^Si state and the 2 ^Sq state has a wavelength of 
1557 nm and an Einstein-A coefficient of 9.1x 10^^ s^^. A focused 1557-nm laser also con- 
stitutes a trap for ultracold atoms in the 2 ^Si state because it is red detuned from the 2 ^Si — > 
2 ^Pj transitions. As the 1557-nm laser light is blue detuned from the 2 ^Sq — t- 2 ^Pi transition, 
atoms in the 2 ^Sq state are anti-trapped. 

experiments Ar^ = 1.01 ±0.13 fm^ {26, 29). Comparing the values of Ar^, we find our result 
to be in good agreement but an order of magnitude more precise. An independent spectroscopic 
measurement in helium on the 2 ^Si — )■ 2 ^Pq transition (30) gives Ar^ = 1.059(3) fm^, ob- 
tained using the most recent QED calculations (19). Although the measurement precision of the 
isotope shift for this transition is comparable to our precision, the smaller uncertainty in Ar^ is 
due to a larger sensitivity to differential nuclear charge effects. Presently, the accuracy to which 
the ^He charge radius is known sets a lower limit on the uncertainty of the '^He charge radius 
determined from helium spectroscopy. Our measurement presents a value for the ^He nuclear 
charge radius of 1.961(4) fm. 

We have also demonstrated that all of the trapped atoms can be transferred to the 2 ^Sq state, 
producing a source of ultracold singlet helium. Optically trapping these atoms simultaneously 
with cold 1 ^So ground state atoms (produced after two-photon decay) opens up the possibility 
to perform two-photon spectroscopy on the 2 ^Sq ^ 1 ^Sq transition {11, 31), where QED and 
nuclear size effects are strongest. 
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Figure 2: Experimental setup. A small fraction of the 1557-nm laser light is split off and cou- 
pled via a fiber-optic link to be referenced to a fiber-based frequency comb. A heterodyne signal 
is monitored on a fast photodiode (PD) to determine the absolute frequency of the 1557-nm 
laser. The remaining light is divided into the trap beam and the spectroscopy beam. A crossed- 
beam dipole trap configuration is realized by focusing both the incident and returning trap beam 
(with orthogonal linear polarizations) to a waist of ~85 /im at the center of the magnetic trap 
(represented by the green coils) under a relative angle of 19 degrees, trapping atoms at the 
intersection. The spectroscopy beam is frequency shifted by a 40-MHz acousto-optical mod- 
ulator (AOM), overlapped with the returning trap beam and absorbed by a thermopile power 
meter (PM). A microchannel plate detector is positioned underneath the trap for temperature 
and atom number determination. 
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Figure 3: (A) Bimodal time-of-flight distribution observed when He* atoms are detected on the 
MCP detector approximately 186 ms after the trapping laser light is turned off. The MCP signal 
is fit to determine the number of Bose-condensed atoms. (B) The percentage of Bose-condensed 
^He atoms remaining in the 2 '^Si state as a function of applied laser frequency (relative to the 
fitted center frequency /o). The line is a fit of a Gaussian to the data. We measure linewidths 
varying from 75 to 130 kHz depending on the trap depth and on the isotope. 
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Figure 4: Measured transition frequencies for ^He (A) and for ^He (B). The error bar on each 
data point includes contributions from various systematics, e.g. AC Stark shift and Zeeman 
shift (23). The frequencies quoted on the y-axes are the weighted averages of the data points. 
Their associated uncertainties (with the uncertainties due to the frequency comb calibration 
and the mean field shift added in quadrature) are ±1.8 kHz for "^He and ±1.5 kHz for ^He, 
represented by the green bands. 
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SUPPORTING ONLINE MATERIAL 



Experimental details 



The 1557-nm laser beam is generated by an NP-Photonics erbium-doped CW fiber laser (Scor- 
pio), which has a specified short-term spectral linewidth of < 5 kHz and has a piezoelectric input 
for fast frequency tuning. The power used for the optical trap is variable and is measured on a 
thermopile power meter. The minimum power required to trap atoms is 15 mW and the maxi- 
mum used is 250 mW, which corresponds to a trap depth ranging from 0.3 to 5 fiK, and axial and 
radial trap frequencies ranging from (cOax, ^rad) = 27r x (20, 140) Hz to 2tx x (80, 570) Hz. The 
lifetime of atoms in the trap is limited by background collisions to approximately 15 s. After be- 
ing loaded into the optical trap, the atoms are illuminated by the spectroscopy beam. This beam 
has a waist of 280 /xm at the position of the atoms, much larger than the cloud size. A measure 
for the number of atoms is obtained by letting the atoms fall onto a microchannel plate (MCP) 
detector, located 17 cm below the trap center. The MCP signal (Fig. 3(a)) shows an integrated 
time-of-flight trace of the atoms as they hit the detector. In the case of ^He the signal shows 
a bimodal character, which is fitted with a linear combination of a gaussian and a parabolic 
function, representing the thermal atoms and Bose-condensed atoms respectively {15, 22). For 
^He, a Fermi-Dirac time-of-flight function is used to fit the MCP signal (77). 



Systematic shifts and uncertainties 



To obtain the value of the transition frequency several systematic shifts need to be evaluated, 
i.e., the AC Stark shift, the Zeeman shift, the recoil shift and, in case of a Bose-Einstein con- 
densate a mean-field shift. A summary of the systematic uncertainties are shown in Table SI. 
The systematic shifts and uncertainties are discussed in detail below. 



Table SI: Combined systematical and statistical uncertainties. Units are kHz. 





^He 


^He 


Frequency comb 


0.4 


0.4 


Zeeman shift 


0.5 


0.5 


Mean-field shift 


1.1 




Recoil shift 


0.0 


0.0 


AC Stark shift extrapolation 


1.4 


1.4 


Total (in quadrature) 


1.8 


1.5 



Frequency comb. The 1557-nm laser frequency is locked to a frequency comb based on a mode- 
locked erbium-doped fiber laser (Menlo Systems), for which both the repetition rate and the 
carrier-envelope-offset frequency are referenced to a GPS -controlled Rubidium clock (PRSIO 
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Stanford Research). The frequency comb has a long-term frequency stabUity of better than one 
part in 10^^. The beat frequency between the 1557-nm laser and one of the comb modes is 
measured on a fast photo-diode. The n*^ comb mode producing the lowest frequency beat with 
the 1557-nm laser is the one for which n-times the repetition rate (frep ~ 250 MHz) of the 
comb laser plus the carrier-envelope-offset frequency (fceo — 20 MHz) matches the 1557-nm 
laser frequency most closely, then /1557 nm — fceo + 'nfrep ± heat- In order to determine the 
absolute frequency of the spectroscopy laser the specific comb-mode number n is determined 
using an accurately calibrated wave meter. The beat frequency of the 1557-nm laser with the 
frequency comb is continuously measured and regulated by a feed-back loop to match a set- 
point frequency. The integration time is 30 ms, which is a trade-off between the averaging 
time of the comb laser and the short-term stabiUty of the spectroscopy laser. This laser locking 
scheme results in a long-term laser linewidth of 75 kHz (FWHM). For each instance of mea- 
suring the 2 ^Si — )■ 2 ^So transition frequency 40 different values of fi,(,at are chosen, where the 
interaction time between the spectroscopy laser and the atoms is, on average, 3 seconds. This 
corresponds to a total interaction time of 120 seconds which results in an expected uncertainty 
of the caUbration with the frequency comb of less than 0.4 kHz. 

Zeeman shift. The transition frequency is shifted by the differential Zeeman shift of the initial 
and the final state in a residual magnetic field, the magnitude of which is approximately 0.5 G, 
predominantly due to Earth's magnetic field. In a field this size, only first-order Zeeman effects 
are of relevance at the precision of our experiment. ^He atoms are initially in the 2 ^Si (mj=-i-l) 
state, introducing a shift rate of gs^sj h = 2.8025 MHz/G, where h is the Planck constant, gs is 
the electron g-factor and is the Bohr magneton. The final state is 2 ^Sq (mj=0), which does 
not shift due to the absence of a magnetic moment. In the case of ^He, both the initial and final 
states have an additional nuclear spin contribution. The initial 2 ^Si, F=3/2 (mF=+3/2) state 
shifts by {gspiB — \gil^N)/h, where gi is the nuclear g-factor and /xat is the nuclear magneton. 
The final state, 2 ^Sq, F=l/2 (mF=+l/2), shifts by —^giUN/h. As such, the differential shift 
rate is gslis/ h, equal to that of ''He. 

For both ''He and ^He, the Zeeman shift is determined by inducing radio-frequency (RF) 
transitions to the mj=0 and mj=-l magnetic substates in ^He (Fig. SI). After a 40-/xs RF pulse 
at a set frequency, the atoms are released from the dipole trap and the spin-state populations 
are separated during time-of-flight in an applied magnetic field gradient. Subsequently, an 
absorption image is taken from which the relative populations of the magnetic substates are 
obtained. The Zeeman shift equals the frequency for which most atoms are transferred. On 
average, the uncertainty of the value we deduce for the Zeeman shift using this procedure is 
0.5 kHz. 

Zeeman shift measurements are performed intermittently with 2 ^Si — )■ 2 ^Sq transition 
measurements in order to compensate for slow changes in the magnetic field. The standard 
deviation of the value of the Zeeman shift over the course of a day is, on average, 2 kHz which 
corresponds to a stability of the magnetic field of ±0.7 mO. 

Mean-field shift. An additional correction to the transition frequency for ^He is the cold col- 
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Figure SI: Measurement of the Zeeman shift. The images show the population of magnetic 
substates of the 2 ^Si manifold after a Stern-Gerlach type experiment as a function of applied 
microwave frequency. The ratio of atoms in the mj=0 and mj=-\ state versus atoms in the 
mj=l state is plotted in the graph. A fit provides the resonance frequency which equals the 
Zeeman shift. 



lision frequency shift (24). Originating from the difference in mean-field energy between the 
ground state and the excited state, this frequency shift is proportional to the density of trapped 
atoms and to the difference in s-wave scattering lengths for the triplet-triplet and the triplet- 
singlet molecular potentials. The triplet-triplet scattering length is -1-7.5 nm (52); however, the 
value of the triplet-singlet scattering length is presently unknown. By measuring the 2 ^Si — )• 
2 ^So transition frequency with a high atomic density and with a low atomic density (reduced 
by a factor of 5) and taking the difference, we find a mean-field shift of 0.07 ± 1.08 kHz which 
is consistent with no mean-field shift to within 1 . 1 kHz. 

Recoil shift. Due to momentum conservation, when an atom is excited from the 2 ^Si state 
to the 2 ^So state, the momentum of the atom is increased by the momentum of the absorbed 
photon: Ap = hf /c, where h is the Planck constant, / is the photon frequency and c is the 
speed of light. The corresponding increase in kinetic energy for the atom must come from the 

photon, leading to a recoil shift of AE = 2^ , where m is the atomic mass. For a ^He 
atom excited by a 1557-nm photon, the recoil shift is 20.6 kHz and for ^He it is 27.3 kHz, with 
negligible uncertainty. 

AC Stark shift extrapolation. The AC Stark shift is linearly proportional to the sum of the power 
in the trap and spectroscopy beams. The trap-to-spectroscopy power ratio is fixed, but the total 
power is variable and measured on a thermopile power meter. The atomic transition is measured 
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Figure S2: Measured transition frequency for '^He as a function of total laser power. The reso- 
nance frequency is AC-Stark shifted as power increases. The field-free transition frequency is 
found at the point where the fitted line crosses the vertical axis. Each data point is obtained by 
performing several measurements and corrected for the Zeeman shift and for the recoil shift. 

several times at a single power and then repeated for a range of powers, as shown in Fig. S2 
for ^He. From a linear fit, the extrapolation to a field-free resonance frequency is deduced. 
On average, one extrapolation gives an uncertainty of 5 kHz. In total, 10 extrapolations are 
determined for ^He and for '^He, reducing the final uncertainties to 1.3 kHz. 

The thermopile power meter readings are corrected for non-linearity by applying a bench- 
top calibration procedure. The uncertainty in this calibration procedure contributes an additional 
uncertainty of 0.6 kHz, leading to a combined uncertainty in the AC Stark shift extrapolation of 
1.4 kHz. 

Theoretical determination of ionization energies 

Indirect determination. Combining experimental transition frequencies from both metastable 
states to higher S, P and D states with theoretical values for the ionization energies of these 
states yields ionization energies 960 332 041.01(15) MHz for the 2 ^Sq state (S, 12, 13) and 
1 152 842 742.97(6) MHz for the 2 ^Si state (2, 8). Combining these values a 2 ^Si ^ 2 ^Sq 
transition frequency of 192 510 701.96(16) MHz is deduced for ^He. 

QED theory. QED calculations of ionization energies of the two metastable states of ^He 
presently agree between different theory groups (S, 9), although their calculations are not fully 
independent. Drake and Yan (S) find 960 332 037.9(1.7) MHz and 1 152 842 741.3(2.5) MHz 
for the ionization energies of 2 ^Sq and 2 "^Si respectively, whereas Yerokhin and Pachucki (9) 
find 960 332 038.0(1.9) MHz and 1 152 842 741.4(2.6) MHz. This implies a theoretical tran- 
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sition frequency of 192 510 703.4(2.6) MHz (for both theory groups), where we took the max- 
imum of the theoretical uncertainties of both ionization energies as the uncertainty in the tran- 
sition frequency since both uncertainties are expected to be correlated (private communication, 
K. Pachucki and V. Yerokhin, 2010). The uncertainty is due to approximations made in calcu- 
lating all higher-order aJ terms in QED, where a is the fine-structure constant. 

Isotope shift. From Table 2 of Ref. {18), the differences in ionization energies between ^He 
and ^He are 45 862 986.0(5) kHz for the 2 ^Sq state and 53 897 130.1(6) kHz for the 2 ^Si 
state. Recently, Drake corrected a sign error in this calculation leading to corrected values 
of 45 862 980.9(7) kHz and 53 897 129.4(8) kHz respectively (private communication Drake 
2010). These values agree with the isotope shifts of the ionization energies as calculated by 
Pachucki and Yerokhin (private communication, 2010): 45 862 979.58 kHz and 53 897 128.33 
kHz respectively. The latter values were calculated applying the methodology of Ref. (9). All 
these values exclude hyperfine structure effects and contributions due to finite nuclear size ef- 
fects, i.e. assuming point charge radii. Taking the difference between the isotope shifts of the 
2 ^So and 2 ^Si ionization energies, the calculated isotope shift for the atomic transition by 
Drake is 8 034 148.5(7) kHz and by Yerokhin and Pachucki is 8 034 148.75(69) kHz, where 
the uncertainty is estimated to be about 20% of the E61 term which is a sum of pure recoil 
and radiative recoil corrections of 0(a^m'^/M), where m is the electron mass and M is the 
nuclear mass. For our calculation of the nuclear shift we use the average of the theoretical iso- 
tope shifts: 8 034 148.63(70) kHz. The hyperfine structure of =^He (25) shifts the 2 ^Si F=3/2 
state by -2 246 587.3 kHz and the 2 ^Sq F=1/2 state by 60.6 kHz (due to hyperfine-induced 
singlet-triplet mixing). 

Helium charge radii from nuclear theory 

Equation 29 from Ref. (18) relates the nuclear charge radius to an effective rms radius of the 
nucleus (corresponding to the distribution of point-like protons and neutrons in the nucleus), 

to the mean-charge radius of the proton and to the mean-charge radius of the neutron. The 
effective rms radius is taken from Tables 7 and 8 of Ref. (28). The results given in these Tables 
correspond to state-of-the-art calculations of nucleon-nucleon and three-nucleon forces. 
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